Quasi-classical cyclotron resonance of Dirac fermions in highly doped graphene 
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Cyclotron resonance in highly doped graphene has been explored using infrared magneto- 
transmission. Contrary to previous work, which only focused on the magneto-optical properties 
of graphene in the quantum regime, here we study the quasi-classical response of this system. We 
show that it has a character of classical cyclotron resonance, with an energy which is linear in the 
applied magnetic field and with an effective cyclotron mass defined by the position of the Fermi 
level m — Ep/vp. 



PACS numbers: 71.55.Gs, 76.40.+b, 71.70.Di, 78.20.Ls 
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I. INTRODUCTION 

The cyclotron motion of charge carriers with a weU de- 
fined effective mass m and the related resonant absorp- 
tion of light at the cyclotron frequency Wc = \e\B/m are 
basic physical phenomena representative of the magneto- 
optical response of conventional condensed-matter sys- 
tems. Well-established and often applied, cyclotron res- 
onance (CR) became through the years a common tech- 
nique routinely used in solid-state physics, inherently 
connected with the effective mass of elementary elec- 
tronic excitations^ 

Recently fabricated graphene^'^ immediately became a 
challenge for CR measurements and implied necessity to 
comprehend the phenomenon of cyclotron resonance in 
a system with massless charge carriers.--'^ Special char- 
acter of CR in graphene can be exemplified on the basis 
of a simple quasi-classical approach.-S"— If we solve the 
classical equation of motion for a charged particle with 
an energy e which depends linearly on the momentum 
p, i.e. e = ±iii?|p|, we find out that the precession-like 
(cyclotron) motion remains characteristic also for these 
Dirac-like particles, nevertheless, with a cyclotron fre- 
quency modified to Wc = |e|i3/(|e|/ti|n). Hence, the classi- 
cal cyclotron motion and consequently also the cyclotron 
resonance of massless Dirac fermions scale linearly with 
the applied magnetic field B, similar to conventional 
massive particles, but they depend on the particle en- 
ergy, which appears in a form of an energy-dependent cy- 
clotron mass defined by the Einstein relation |e| = mvp. 

Surprisingly, this classical, i.e. linear with the mag- 
netic field, cyclotron resonance has not been clearly ob- 
served in graphene so far. Instead, the optical response 
of graphene has been experimentally explorediSrJ^ only 
in the quantum regime, when all detected transitions cor- 
respond to excitations between well-separated Landau 
levels (LLs). These measurements provided us with a 
unique insight into the LL spectrum of graphene and di- 



rectly visualized its intriguing yS-scaling, but no infor- 
mation about the quasi-classical limit could be extracted. 
Consequently, it remains a challenge to demonstrate ex- 
perimentally the quasi-classical cyclotron resonance in 
graphene, or interestingly, to study the crossover between 
this classical (linear in B) and fully quantum-mechanical 
(linear in VB) magneto-optical response of graphene - 
the crossover, which is not present in a system of mas- 
sive particles, where the response always remains linear 
in the applied magnetic field. 

Here we report on magneto-optical measurements on 
highly-doped epitaxial graphene grown on the silicon- 
terminated surface of SiC. We identify the cyclotron- 
resonance-like response related to the massless Dirac 
fermions in the vicinity of the Fermi level, which evolves 
nearly linearly with the applied magnetic field, as ex- 
pected from (quasi-)classical considerations. These re- 
sults are compared to the currently well-understood re- 
sponse of multilayer epitaxial graphene prepared on the 
carbon-terminated surface of SiC, where majority of (ro- 
tationally stacked) graphene sheets is nearly undoped 
and where the fully quantum-mechanical treatment is ap- 
propriate for its description. 



II. SAMPLE PREPARATION AND 
EXPERIMENTAL DETAILS 

In our study, we present data taken on two particular 
samples of graphene prepared by thermal decomposition 
on the surface of semi-insulating silicon carbide^ Both 
specimen have been grown at 1600°C in Ar atmosphere 
using a commercially available horizontal chemical vapor 
deposition hot-wall reactor (Epigress V508), inductively 
heated by a RF generator. 

The sample denoted as A in our study, grown on the C- 
terminated surface of silicon carbide (4II-SiC[000i]), rep- 
resents a standard multilayer epitaxial graphene (MEG) 
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FIG. 1: (color online) Relative transmission spectra of the 
sample A (MEG on C-face of SiC) for selected values of the 
magnetic field. For clarity, successive spectra are shifted ver- 
tically by 0.1. Individual transitions discussed in this paper 
are denoted by vertical arrows. 



specimen. The majority of sheets in this structure is 
practically undoped (residual p-like doping expected at 
the level of 10^ cm~^, see Refs. |l3| and 1 141) with the elec- 
tronic band structure similar if not identical to an iso- 
lated graphene monolayer, provided by the characteristic 
rotational stacking of layers. ^^ This fact is reflected in 
the micro-Raman spectra taken on this sample, which 
revealed dominantly single-component (Lorentzian-like) 
2D bandji?'^^ nevertheless, some Bernal-stacked residu- 
als have been observed on selected locations as well. 

The sample B has been prepared on the Si-terminated 
surface of SiC (4H-SiC[0001]) and also characterized in 
the micro-Raman experiment. Similarly to the previ- 
ous sample, the dominantly observed single-component 
2D band serves as an indication for linearity of elec- 
tronic bands. It points towards the presence of a sin- 
gle graphene sheet (in addition to the non-graphene-like 
buffer layer—), as in contrast to MEG on the C-face, the 
few-layer graphene stacks on the Si-face exhibit Bernal 
and not rotational ordering.— In addition, such mono- 
layers are typically rather highly doped by the charge 
transfer from the substrate i^i^i 

To measure the (non-polarized) infrared transmittance 
of our samples, we used the radiation of a globar, which 
was analyzed by a Fourier transform spectrometer and 
delivered to the sample via light-pipe optics. The light 
was then detected by a composite bolometer kept at 
T = 4.2 K placed directly below the sample. Measure- 
ments were done in the Faraday configuration in a su- 
perconducting solenoid. Both investigated samples were 
fully opaque in the energy range 85-220 meV due to 
phonon- related absorption in SiC. 




4 6 8 

Magnetic field (T) 



10 



FIG. 2: (color online) Positions of absorption lines observed 
in the A sample grown on C-terminated surface of SiC. The 
spheres correspond to inter-LL transitions in electrically iso- 
lated and nearly undoped layers in MEG, having a charac- 
teristic v'^-dependence. The dashed lines show their theo- 
retically expected positions for the Fermi velocity of vp = 
1.02 X 10® m.s~^. Origin of the transition denoted by stars, 
which evolves nearly linearly with B, is discussed in the text. 



III. RESULTS AND DISCUSSION 

Let us start our discussion with the data obtained on 
the sample A, i.e. on the MEG structure grown on the 
C-terminated surface of SiC, see Fig. [U whose magneto- 
optical response is relatively well understood] ^°'^'^i^^ 
Transmission spectra are in this kind of specimens dom- 
inated by relatively sharp absorption lines with a typ- 
ical \AB-hke scaling (see Fig. [5]) and taking the LL 
spectrum of graphene: En — s\gn{n)vp \/2h\ eBn\, we 
easily identify them as inter-LL transitions L_„— >L„+i 
and L_(„_|.i)-^-L„ with n — 0...6. The corresponding 
Fermi velocity can be evaluated as vp = (1.02 ± 0.01) x 
10^ m.s~^. As a matter of fact, this response is equiv- 
alent to that of an exfoliated graphene monolaye r^^'^^i^^ 
and in turn it serves as a direct indication that layers in 
MEG indeed behave as independent, electrically isolated, 
graphene sheets. 

Apparently, here we deal with graphene in its fully 
quantum regime, i.e. we study optical excitations be- 
tween well-separated Landau levels, whose -s/B-scaling 
and multi-mode character are directly related to the lin- 
earity of the electronic bands*?'* To characterize our sam- 
ple further, we can set the upper limit for the carrier den- 
sity that cannot significantly exceed 10^° cm~^, as the 
transition Lo(-i) ~^L]^(o') is observed down to the mag- 
netic field of 100 mT (the filling factor v < 6). We as- 
sume that the graphene sheets are slightly p— doped as 



follows from the recent ARPES study done on equiva- 
lent specimens J^ We should also note that the dominant 
graphene-like lines are often accompanied by additional 
spectral features, which are usually significantly weaker, 
nevertheless still well-defined in the spectra, located es- 
pecially at low energies, see e.g. the transition denoted 
by stars in Fig. [TJ We will come back to these features 
later on after we discuss the response of the sample B. 

The sample B exhibits a very much different behav- 
ior in our magneto-transmission experiment and shows 
no vS-scaled spectral features, which are, as we have 
learned above, characteristic of graphene in its quantum 
regime. Instead, it clearly shows magnetic-field-induced 
transmission at low energies accompanied by a relatively 
wide absorption minimum, which moves towards higher 
energies with increasing B. We found this behavior to 
be characteristic for all graphene specimens prepared on 
Si-terminated surface of SiC (not shown in this paper), 
nevertheless, the total intensity of the observed spectral 
features as well as the position of the absorption mini- 
mum slightly varied from sample to sample. 

To explain this behavior, we should recall that epi- 
taxial graphene on the Si-terminated surface is usually 
highly doped by the charge transfer from the substrate, 
reaching carrier densities up to 10^^ cm^^ {Ep w 0.4 eV), 
see e.g. Ref. [20, with quality (expressed in terms of scat- 
tering time or mobility) comparable to MEG on the C- 
face of SiCJ^i^^ Therefore, at magnetic fields of a few 
Tesla, LLs should not be resolved around the Fermi 
level and the system should remain in the quasi-classical 
regime. Interestingly, owing to the non-equidistant spac- 
ing of LLs in graphene, the tuning from a fully quantum 
regime to the classical one is possible just by changing 
the Fermi level. 

If we adopt this, the qualitative explanation of the ob- 
served behavior becomes straightforward. The relative 
magneto-transmission spectrum T{B)/T{0) reflects two 
specific effects: The low-energy free carrier (Drude-like) 
absorption, which is gradually suppressed with the in- 
creasing field, and the cyclotron resonance itself. The 
former one gives rise to the field-induced transmission at 
low energies that decreases monotonically with the pho- 
ton energy, and the latter one is manifested by an ab- 
sorption minimum close to the energy of the cyclotron 
resonance ujc- Let us now proceed with the quantitative 
analysis of this data, which will allow us to determine the 
characteristic field- dependence of Wc. 

In the framework of the classical theory for cyclotron 
resonance, the optical conductivity is given by 
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which describes both the CR active (— ) and inactive (-I-) 
polarization modes of light, as well as the Drude-like low- 
energy absorption (at -6 = 0). ctq here denotes zero- field 
dc conductivity and 7 is a phenomenologically introduced 
damping. We should also point out that in graphene, 
contrary to conventional systems of massive particles, the 
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FIG. 3: (color online) Transmission spectra of the sample 
B taken at various magnetic fields in the far infrared spec- 
tral region. Sohd curves represent fits to experimental data 
(spheres) according to the quasi-classical theory described in 
the text. For clarity, successive spectra are shifted vertically 
by 0.1. Vertical arrows indicate position of a relatively broad 
CR resonance, as obtained from the analysis. The inset shows 
fitting procedure at -B = 10 T in detail. The individual 
transmission curves T{Q),T{B),T+{B) and T-{B) have been 
plotted to demonstrate their contribution to the final relative 
change of sample transmission T{B)/T{Q). Let us note that 
these curves share a common baseline given by the transmis- 
sion of the slab without any conducting layer (or away from 
the resonances) Tsiab = 2n/(n^ -|- 1) « 0.67. 



cyclotron frequency lUc is not independent of ctq a-nd 7 as 
the "effective" effective mass of Dirac fermions (at the 
Fermi level) can be determined from these two parame- 
ters. 

The transmission of circularly-polarized light through 
the studied specimen (graphene sheet on an isolating sub- 
strate) can be expressed by formula)^ 



T±{io,B) = 
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which is derived for a dielectric slab with the index of 
refraction n {nsic — 2.6) with an infinitely thin con- 
ducting layer on one of its surfaces (described by cr± ) . c 




2 4 6 8 10 12 

Magnetic field (T) 



FIG. 4: (color online) Position of the cyclotron resonance 
in the sample B as extracted from the fitting experimen- 
tal curves in Fig. |3l well-following the linear in B depen- 
dence plotted for the cyclotron mass of m, = O.OSSttio (full 
line). For comparison, we also depict positions of intraband 
(i.e. cyclotron) resonances between LLs with higher indices, 
L„ — >-L„+i n = 6 ... 40 (doted lines). The inset schematically 
shows linear in B evolution of CR in highly doped graphene, 
despite the fact that the individual LLs exactly follow the 
vB-dependence. As we cannot determine the sign of charge 
of the particles in our non-polarized transmission experiment, 
we considered n-type doping only in this inset for simplicity. 



and £o stand for the real speed of light and the vacuum 
permittivity, respectively. The dependence of T± on the 
thickness of the slab (i.e. interference effects) has been 
aheady averaged out. 

Within the framework of our simple model, we can 
easily calculate the relative change of the sample trans- 
mission with the magnetic field [T+{B) + T_(B)]/2T(0), 
and compare it with our experimental data. As the con- 
ductivity as well as the damping parameter can in prin- 
ciple vary with B, we have first found their appropriate 
values at zero field, ao = (2.4 ± 0.3) x 10~^ n~^ and 
7o = (10 ± 2) meV, which we afterwards kept constant 
and varied only Wc, crs and 73 to reproduce our experi- 
mental traces. 

As shown in Figs.[3]and|ll our simple theoretical model 
is able to reproduce experimental transmission spectra 
very well and the CR frequency indeed follows a clear 
linear in B dependence. The conductivity ctb and the 
CR half-width jb exhibited in our fitting procedure a 
weak decrease and increase, respectively, but practically 
remain within the interval of uncertainty estimated for ctq 
and 7o and they can be therefore taken as nearly constant 
in the investigated range of fields. The width (FWHM) 



of the cyclotron resonance absorption peak 27^ has been 
found comparable with the cyclotron energy Wc what jus- 
tifies the quasi-classical regime in the investigated sam- 
ple. 

Wc should also clarify that a relatively broad CR res- 
onance implies that the minimum in the transmission 
spectra does not exactly match to the position of the 
CR mode. This is explained in the inset of Fig. [3l 
where the transmission curves T{B) at _B = 10 T and 
T(0) are depicted as calculated from parameters obtained 
in our data analysis. We also plot curves T-{B) and 
T-\-{B), which show expected transmission in opposite 
circular polarizations of light, i.e. in the CR active 
and inactive modes, respectively. Let us note that all 
these transmission curves share a common baseline at 
2n/{n^ -\-l) K, 0.67, which corresponds to the light trans- 
mission through a dielectric slab without any conducting 
sheet (or away from resonances). This value can be easily 
obtained from Eq. ^ for a = Q. 

Let us now return to the main finding of our anal- 
ysis, i.e. to the clear linear in B dependence of the 
cyclotron resonance as shown in Fig. |4l If we evalu- 
ate the effective mass corresponding to this CR mode, 
we obtain m — (0.053 ± 0.004)mo, which according 
to the Einstein relation Ep — mvp implies the Fermi 
level of Ep « 310 meV in a reasonably good agree- 
ment with conclusions of ARPES measurements on sim- 
ilar samples j ^°'^^ In this estimation, we used the Fermi 
velocity of 1.02 x 10^ m.s~^ determined from the mea- 
surement on the sample A, but very similar velocities are 
reported by ARPES for graphene on the Si-terminated 
surface of SiC as welh ^°i^^ 

The observation of the CR mode that is linear in the 
applied field B justifies our simple approach in which 
(otherwise massless) Dirac fermions are considered as 
massive particles, the mass of which is determined by 
the Fermi energy. On the other hand, this confirma- 
tion implies another constrain on parameters used in 
our fitting. Namely, the conductivity ctq can be di- 
rectly estimated within the Boltzmann transport theory 
from the scattering rate 70 /t and from the Fermi en- 
ergy: CTo — {2e'^/h){Ep/jo). This relation was respected 
in our fitting procedure and required an additional self- 
consistent loop to keep the correct relation between ctq 
and 7o on the basis of Ep extracted from the fitting. 
Obviously, this approach works only in case of specimens 
with a full coverage by the graphene monolayer, other- 
wise the conductivity derived from 70 and Ep exceeds 
the value obtained directly from the fitting. The ratio of 
these two values thus could serve as a simple parameter 
that describes fraction of coverage by graphene. 

We should note that the evaluated mass of 0.0537tt.o 
is surprisingly close the effective mass of K point elec- 
trons in bulk graphite ~0.058morl but this similarity 
represents just a coincidence. To prove this, we can for 
instance argue that bulk graphite reaches the quantum 
limit at fields around 7 T^^^ and therefore it cannot 
exhibit CR line that corresponds to the effective mass 



of K point electrons. Instead, due to the formal sim- 
ilarity of the band structure at the K point with the 
graphene bilayer, it should provide the CR-like absorp- 
tion at energy by a factor of \/2 higher, as indeed ob- 
served experimentally;'^^ 

To summarize, the highly-doped graphene sheet on the 
Si-terminated surface allows us to study the CR reso- 
nance of masssless Dirac fermions in its fully classical 
limit, when it becomes strictly linear with the applied 
magnetic field. Nevertheless, let us look at the same 
problem using the quantum-mechanical description, i.e. 
interpret our data as intra-band excitations between ad- 
jacent (but in reality significantly broadened) Landau 
levels, as schematically shown in the inset of Fig. 21 To 
clarify the linear in B response, found in a system with 
a vi?-scaled LL spectrum, we only need to consider the 
Fermi level Ep located high in the conduction (or equiv- 
alently low in the valence) band and independent of the 
magnetic field. This can be exactly ensured by the pin- 
ning the Fermi level by the surface states in the substrate 
or fulfilled approximatively, when LLs significantly over- 
lap. At a given magnetic field, the n-th LL approaches 
the Fermi level, Ep = Ei^/n, and the active CR-like 
transition (between adjacent LLs) has the energy of 
hw, = Si(V^rTT- V^) « Ei/{2^) = heB/iEp/vl), 
i.e. it indeed becomes linear in B and moves with the 
slope related to the cyclotron mass m — Ep/vp. 

The analysis of data obtained on highly-doped 
graphene in the sample B allows us to come back to the 
low-energy features in spectra taken on the sample A. 
Also in this kind of samples, we usually find relatively 
weak spectral features whose energy moves nearly lin- 
early in B. Even though the sample A dominantly con- 
sists of nearly undoped graphene sheets, the increase of 
the magneto-transmission at lower energies serves as an 
indication of some free carriers in the sample, similarly to 
the sample B. In fact, also in MEG on the C-face of SiC, 
first few layers nearby the interface are usually doped by 
charge transfer from the substrate. The transition de- 
noted by stars (see Figs. [1] and [S]) thus could be assigned 
to the CR response of these doped layers, and its charac- 
teristic slope «4.5 meV.T^^ (effective mass of 0.025too) 
would imply the Fermi level of \Ep\ « 150 meV. On 
the other hand, the carrier density in these doped lay- 
ers gradually decreases with the distance from the sub- 
strate, which should prevent observation of a well-defined 
CR mode (whose energy depends on the density of mass- 
less particles), while the transition denoted by the star 
remains relatively narrow in our data. 

More likely, this line originates in Bernal-stacked resid- 
uals, especially in the graphene bilayer which has been 



identified as a minor but not negligible component 
of MECJii^ The masssive Dirac fermions in the bi- 
layer are characterized by an effective mass of ms = 
li/i^vl) « 0.32TOO (71 = 380 meV^-^) which does 
not match to the effective mass derived for the "star" 
transition, but the characteristic LL spectrum in the 
bilayer, _B„ = sign(n)(/i|ei3|/mB)-y/|n|(|n| -I- l)^iS im- 
plies in the quantum limit appearance of the CR mode 
(L_i(o) ^•Lo(i) transitions) whose slope is by a factor 
of \2, higher. The transition denoted by the star, with 
a slope corresponding to the effective mass of 0.025mo, 
thus could be assigned an undoped graphene bilayer in 
which the CR absorption should correspond to the mass 
of ttib/v^ w 0.023?tt,o. Nevertheless, the definitive con- 
clusion cannot be drawn on the basic of our current data. 



IV. CONCLUSIONS 

We have explored magneto-optical response of 
graphene in two distinctively different regimes. Whereas 
the weakly doped sheets in multilayer epitaxial graphene 
on the C-terminated surface allowed us to study the opti- 
cal response of graphene in the quantum limit, massless 
Dirac fermions in the quasi-classical regime can be ex- 
plored in highly doped epitaxial graphene prepared on 
the Si-terminated surface of SiC. The response of the 
former system clearly shows vB-like scaling, typical of 
Landau levels in systems with linear electronic bands, 
the latter case proves that cyclotron resonance of mass- 
less particles remains in the quasi-classical limit linear 
with the magnetic field. 

Additional note: During our work on the manuscript we 
became aware of Faraday rotation measurements on epi- 
taxial graphene, ^^ in which a similar conclusion about 
cyclotron resonance of massless Dirac fermions in the 
quasi-classical limit has been made. 
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